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ABSTRACT 


Although acoustic detection of cavitation inception has 
been shown to agree relatively well with visual detection, 
acoustic methods have generally not been used to detect 
cavitation inception during cavitation testing of model 
propellers. In addition, it has been suggested that noise 
Measurements on model propellers be made at high frequencies 
to more properly represent the full scale noise. In this 
thesis, three different methods of acoustic detection were 
investigated. Two of these methods, the measurement of high 
frequency one-third octave band levels and the analysis of the 
complete noise spectrum between 10 and 50 kHz, met with some 
success, but were not equivalent to the capability of a 
visual detection method. The third method used, the demodulated 
analysis of high frequency cavitation noise, gave excellent 


agreement with visually determined results. 
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ales INTRODUCTION 


Since the first observation of cavitation associated with 
Marine propellers was reported by Reynolds in 1873, a number 
of unwanted effects, including loss of propeller efficiency, 
erosion of propeller surfaces, excitation of hull vibrations, 
and generation of noise, have been identified and studied. 
Because of these detrimental effects, it has been, and 
continues to be, desirable to predict the cavitation 
performance of a propeller design before the expensive full 
scale propeller is built. With no exact analytical approach 
available for predicting the full scale cavitation performance 
of a propeller, the testing of scale models has been used to 
aid in cavitation prediction in the propeller design process. 

For the model test to properly represent the full scale, 
it 1s necessary for Similarity conditions be satisfied. For 
propeller cavitation testing, this amounts to uSing a 
geometrically similar propeller operating in a flow which 
matches the wake where the full scale propeller operates. 

With these conditions met, it is assumed that cavitation 
performance for similar values of cavitation index, 
yoy ee 


ko US 


ane advance coefficient, 
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will be the same for the model and full scale propeller. 
But this assumption is not precisely correct, and scale 
effects, which arise from the inability to satisfy all 
hydrodynamic, thermodynamic and other microscopic similarity 
requirements, are encountered. These scale effects are 
usually eliminated by means of empirically or theoretically 
determined corrections. 

In general, the procedure for conducting a model test for 
determining cavitation inception performance is to operate 
the model propeller in a variable pressure water tunnel, 
downstream of a device which produces the desired wake at the 
plane of the propeller. A water and propeller speed combination 
are chosen to give the desired value of advance coefficient. 
Water pressure is changed to change the cavitation index 
until cavitation is visually observed to either begin or to 
cease, depending upon the criterion used at the particular 
test facility. This process is then repeated for several 
different values of advance coefficient. The final result is 
a curve of inception cavitation index, Os versus advance 
coefficient. 

Although visual observation is the usual method for 
determining the presence or absence of cavitation, it 1s not 
the only available means. It is possible to use the detection 
of cavitation-generated noise to determine, or to assist 
visually determining, the inception of cavitation. Good 
correlation between acoustic and visual inception determination 


has been reported ( Lehman, 1964). It has also been reported 





-10- 
that "numerous" facilities use an acoustic technigue for 
this purpose (ITTC, 1978), but the details of these methods 
were not available. It is proposed for this investigation 
to consider several different schemes for detecting cavitation- 
generated noise as a method for inception determination, and 
to compare acoustically-determined inception data with visual 


observations of inception for a model propeller. 
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If. BACKGROUND 

The sound generated by cavitation comes primarily from 
the growth and collapse of the cavitation bubbles. The 
theoretical energy spectrum for the sound generated by a 
Single bubble has been shown to contain maxima at frequencies 
which correspond to the reciprocal of the time required for 
the growth and collapse of the bubble (Fitzpatrick and 
Strasberg, 1959). Experimental investigations into the 
spectrum of cavitation noise have found that the shape of the 
measured spectrum resembles the theoretical spectrum quite 
closely (Ross, 1976; Strasberg, 1977). 

Strasberg also notes that the peak of the observed 
spectra move toward lower frequencies as the cavitation 
becomes more severe. The larger maximum size of the bubbles 
in the more developed cavitation corresponds to the observed 
peak at a lower frequency. ROSS points out that the energy 
radiated per collapse is proportional to the product of the 
collapse pressure and the maximum bubble volume. So, when 
cavitation becomes more severe, and a larger number of 
bubbles,which also have a greater diameter, are produced, the 
amount of sound energy radiated becomes greater, the magnitude 
of the peak in the spectrum increases, and the frequency of 
the peak becomes lower. If the various spectra are non- 
dimensionalized in the manner of Fitzpatrick and Strasberg, 
the spectra for different degrees of cavitation intensity 
all agree well with the non-dimensionalized theoretical 


spectrum (Strasberg, 1977). 





= 

Continuing on with the scheme given by Strasberg, if the 
propeller diameter is used for the non-dimensionalizing length 
scale instead of the maximum bubble radius, a similarity 
condition for relating frequencies of interest between a 
model and a full scale propeller is obtained. For a given 
ratio of maximum bubble radius to propeller diameter (which 
can be interpreted as a measure of relative cavitation 
intensity), the non-dimensional frequency of the peak in the 
cavitation noise spectrum will remain invariant between 
different length scales. It should then be possible to 
compare cavitation noise measurements made at a given actual 
Frequency on a full scale propeller with measurements made at 
the same non-dimensional frequency on a model propeller, so 
long as other similarity requirements are satisfied. For 
example, with the submarine propeller cavitation yer 
measurements used by Strasberg (1977), assuming that the full 
scale measurements were made at a submergence depth of 200 
feet, the model measurements were made at an ambient pressure 
of 1 atmosphere, the same fluid was used in each case, and 
the length scale ratio was 8, the ratio between the frequency 
of interest with the model (fi) to the frequency of interest 


with the full scale propeller (f) 1s: 


f p )? 
M Dy M 
fo Duy Po 





L 
8° (34 + 234)7 = 3.05 
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So, for cavitation noise measurements, the frequency of 
interest in model scale is about three times the full scale 
frequency for equivalent severity of cavitation in the two 
cases. 

There are several considerations about the model cavitation 
noise which this concern with frequency alone does not show: 
1) The size of the cavitation bubbles for the model 
Will be one eighth the actual size of the full scale 
bubbles and might, consequently, be too small to see 

2) At the same distance from the propeller as in 
model measurements (assuming the distances are large 
enough to avoid near field effects) the sound vressures, 


D, would have the ratio of 

















N 
Pm MM 
——_ 
Pp uP Pp 
ee 10° 
234 8 


Or approximately 35 db lower for the model, if the 
non-dimensionalization of Strasberg (1977) is used with the 
Same bandwidth, distance and non-dimensional frequency. 
3) The actual cavitation does not, in general, 
occur uniformly for all angular positions for all 
blades. 
Operation of the propeller in a non-uniform wake causes 
variations of inflow velocity seen by the propeller which are 
periodic, with a frequency that corresponds to once per 


revolution. This periodic flow variation causes the inception, 
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growth, decrease, and disappearance of cavitation to occur in 
a periodic fashion on a given blade. For a given average level 
of cavitation intensity, the amplitude of the cavitation noise 
will vary over one propeller revolution. This change in 
amplitude will have two effects - it will shift the peak 
frequency of the noise spectrum over the time span of one 
revolution of the propeller, and it will vary the amplitude 
of the noise spectrum. 
If all blades of the propeller are the same, the 
amplitude modulation of the noise occurs at blade passing 
rate - once for each blade for each revolution of the shaft. 
Since the blades are generally not identical, one blade will 
usually begin to cavitate ahead of the others, and the 
MeGgdueltton Of the noise will occur, in addition, at the shaft 
rate (Strasberg, 1946; Ross, 1976). 
It is intended, then, to investigate the use of these 
two aspects of cavitation noise, frequency scaling and 
amplitude modulation, either independently or together, as a 
means of detecting cavitation inception on a model propeller. 
It is expected that this approach would have certain advantages 
as a part of the process for predicting full scale cavitation 
performance: 
(1) Visual determination of inception is very dependent 
upon a number of conditions outside the test tunnel for 
repeatable results. Lighting conditions, ae well as the 


location and visual acuity of the observer, can have a 
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substantial effect upon the outcome of a test. With an 
appropriate criteria for determining cavitation inception 
from acoustically obtained data, this sort of variation 
could be eliminated. 

(2) For many ships, cavitation inception determination 
for the full scale propeller is accomplished using 
acoustic information. An acoustic method on model scale 
would more closely approximate the full scale test. 

(3) Based upon bubble size considerations, an acoustic 
method might be able to detect the presence of cavitation 
bubbles before they are visible. 

But there are disadvantages associated with using 
acoustic information for this purpose: 

(1) The equipment to make the acoustic measurements 
is substantially more expensive than that needed for the 
visual determination of inception. For a very large 
length scale ratio, the frequency of interest at the 
model scale might become so high that the normal analysis 
equipment for acoustic measurements would not be usable, 
or the level of the acoustic signal from cavitation noise 
might be too low to be detected. 

(2) If it is used alone, an acoustic method would 
appear to be less useful to the designer, since the 
method would not directly identify the type of cavitation 
causing the noise. The steps necessary to improve the 
cavitation performance of an unsatisfactory design would 


be less clear. 
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1B lig lag EXPERIMENTAL PROCEDURE 


ne Equipment setup 


Experiments were conducted in the MIT Variable Pressure 
Water Tunnel, which has a 20 in square closed jet test section 
with a length of 54 in. All four walls of the test section 
have a 16 in by 44.5 in plexiglass viewing window insert. 

The propeller is located at the vertical and horizontal 
centerlines of the test section, and is driven by an upstream 
propeller shaft. 

The propeller used for this test was the David Taylor 
Naval Ship Research and Development Center (DTNSRDC) model 
3927, which had a diameter of 10.8 in and had seven blades. 
The tapered hub of the propeller had a maximum diameter of 
2.8 in, which required the use of an adapter to provide a 
smooth transition between the 2.375 in diameter of the 
propeller shaft housing and the hub. Figure 1 shows this 
adapter. The propeller was installed on the shaft with no 
hub fairwater cap, leaving the mounting capscrew and 


lockwasher exposed. 


Attached to the shaft housing 20 in upstream of the 
plane of the propeller was the holder for the screen which 
generated the desired wake at the propeller. The design of 
this screen used a scheme proposed by McCarthy (1963) as a 
Starting point. The details of this process are given in 
Appendix A. 


The wake prescribed for testing this propeller is 
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axisymmetric, with a specific radial distribution of 
longitudinal velocity. For the first variation of the wake, 
the velocity distribution was measured along a diameter on a 
diagonal, but at the radii where the values of longitudinal 
velocity were given. The values for the two radii were 
averaged and compared to the specified values. An error of 
less than 10-12% was considered acceptable. 

During the initial testing of the propeller it was 
discovered that face cavitation would occur behind the 
supports of the wake screen holder. This indicated a velocity 
increase as the propeller blade entered the region downstream 
of the supports, and was attributed to boundary layer viscous 
effects as the flow passed the wake screen holder supports. 
At the same time it was noted that modulation of the 
cavitation noise for other than this face cavitation was not 
detectable with the equipment being used. This indicated 
that a severe, once per revolution, velocity defect was 
desirable. This defect was achieved by using screen material 
to make the topmost support much thicker and tapered. No 
effort was made to maintain the same circumferential mean wake. 
The final wake screen used is shown in the photographs in 
figures 2A and 2B. Figure 3 shows the results of the wake 
survey made with a 1 in square grid, in a plane 2.5 in 
downstream of the blade root leading edge, but with the 
propeller removed, after the upper support was altered. The 
diagonal line indicates where the initial screen velocity 


measurements were made, as well as the data from the final 
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Final Wake Screen - Looking Downstream 


Figure 2A 
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Final Wake Screen - Looking Upstream 


Figure 2B 
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screen velocity measurements, for comparison with the 
prescribed wake were taken. A plot, comparing the actual and 
prescribed velocities is shown in figure 4. 

Two sensors were used at different times to obtain the 
acoustic signal. The first, an accelerometer, was a Briel & 
Kjaer (B&K) type 4344. The characteristics of this 
accelerometer are shown in table 1. The accelerometer was 
mounted directly to one of the viewing windows, as close to 
the center as possible, using a cyanoacrylate adhesive. 

The other sensor, a Minature hydrophone, was a B&K type 3103. 
The principal characteristics of this hydrophone are given in 
table 2. The hydrophone was mounted in the viewing window 

as shown in figure 5, 2.5 in downstream of the leading edges 
of the propeller blades. A schematic diagram of the 
arrangement of the test section is shown in figure 6. 

The methods of processing of the signal from the sensor 
are shown in figures 7A and 7B. In configuration 7A, the 
Ithaco 4213 filter was used as a hand pass filter for one- 
third octave bands, with the level indicated on the B&K type 
2607 measuring amplifier as the output. The other 
configuration used a Federal Scientific Model UA-15A 
Ubiquitous Spectrum Analyzer coupled to a model 1015 
Spectrum Averager, which, in turn, drove an X-Y¥ plotter 
to provide an output. 

With the spectrum analyzer providing the output, two 
set-ups were used. The first was to obtain the complete 


spectrum of the cavitation noise to 50 kHz, the upper 
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Wake Fraction versus Non-Dimensional Radius 


Figure 4 
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Table l 


Accelerometer Characteristics 


Type: Bé&K 4344 Serial Number: 475507 
Reference sensitivity at 50 Hz at 23 °c and including cable 
capacitance of 106 pF: 


Voltage sensitivity: 0.308 mV/ms~* 


Gf o 02am /¢ 

Charge sensitivity: 0.344 pC/ms ° one Sree jolS yee 
Capacitance (including cable): 1116 pF 
Weight: 22 gm Undamped natural frequency: 121 KHz 


Frequency response: 


30 
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Table 2 
Hydrophone Characteristics 
Type: B&K 8103 Serial Number: 636764 


Reference sensitivity at 250 Hz at 23 2 including 6m 


integral cable: 


Open circuit sensitivity: 


Voltage sensitivity: -211.6 db re 1V/uPa 
=91.6 Gb te 2V/Pa or 26.3 lV per 
Pa 
-1l11.6 db re 1V per ubar 
Charge sensitivity: Dil SSG —— pC per Pa 


Capacitance: 3480 pF 


Frequency response: 
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Figure 6 
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frequency limit of the spectrum analyzer. In this case, the 
4213 filter was used as a 10 kHz high pass filter to prevent 
the large amount of noise below 10 kHz, little of which was 
considered to be cavitation related, from overdriving the 
analysis equipment and preventing a good representation of the 
high frequency noise. For this analysis, 128 spectra were 
averaged to provide the output. 

The second setup was intended to detect the cavitation 
by the presence of modulation of the high frequency noise, so 
that some form of demodulation was required. Demodulation of 
the signal was performed by passing it first through the 
Bolt, Beranek and Newman (BBN) Transient Signal Analyzer, 
which squared the signal, and then through a 5 kHz low pass 
. Filter and into the spectrum analyzer. For this setup, the 
4213 filter was used either aS a 20 kHz high vass filter or a 
50-63 kHz band pass filter, so that only the high frequency, 
Ccavitation-related noise was being analyzed. (It was 
necessary to use high frequencies to obtain good results with 
this method of demodulation.) The spectrum analyzer oe used 
on the 0-500 Hz range, and 32 spectra were averaged to obtain 


Ene OUtLPUt. 
Peeeeecalabration 


Tt was not intended to attempt to measure the absolute 
levels of the cavitation noise, so a calibration of the level 
of the signal was not performed. In addition, a calibration 


of the frequency display of the spectrum analyzer was not 
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performed. Consequently, many of the spectra of the 
demodulated signal do not have peaks at the proper 
frequencies. But, because the peaks generally showed the 
proper spacing, this error in the display was not considered 


Sime faucant . 


Cc. Test Procedure 


The conventional test procedure of maintaining constant 
free stream velocity and propeller speed, and varying static 
pressure was used for this testing. This technique kept the 
frequency of interest for the demodulated analysis constant 
for a given test run. 

Since it was not possible to determine the air content 
of the tunnel water, it was decided to perform the acoustic 
and visual cavitation inception determinations concurrently. 
Thus variations in air content would affect visual and 
acoustic results equally. 

Prior to beginning data recording for a series of data 
points, the water was drained from the test section and the 
propeller was operated in air at about 1000 RPM. At this 
time, tare loadings for the thrust and torque load cells were 
recorded and the height of the manometer column was recorded 
as a no flow condition zero. The test section was refilled 
with water, an..initial atmospheric pressure reading was 
recorded, and testing began. 

For each data point, the following sequence was 


followed: 
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(1) A nominal flow speed and propeller RPM were 
selected to give the desired value of J, and tunnel conditions 
were adjusted to these values. The flow speed chosen was such 
that cavitation inception would occur after tunnel static 
pressure was reduced from atmospheric pressure, but before 
the SCORE was so low that air coming out of solution would 
begin to cause absorption of the high frequency acoustic 
Signal or to obscure the propeller from view (250-300 mm Hg). 

(2) Tunnel static pressure was set to atmospheric 
pressure. Room air and tunnel water temperatures, and 
amplifier gain and filter settings were recorded, and the first 
set of data taken. This data included tunnel static pressure, 
manometer height, propeller RPM, thrust and torgue readings, 
and the acoustic data, either the one-third octave level or 
the spectral analysis. A visual observation of the propeller 
was made uSing a strobe light triggered by a once per 
revolution signal from the propeller shaft. A variable 
triggering delay on the light allowed the observation of all 
blades at all points in the propeller rotation. The large 
viewing windows enabled viewing the propeller under 
conditions of both back and front lighting, from both up and 
downstream using the one strobe light. 

(3) Tunnel static pressure was lowered, and another 
set of data recorded. This process was continued until 
cavitation inception had been observed both visually and 
acoustically. 


The sequence above was then repeated at the selected 
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values of advance coefficient needed to produce the curve of 
Cavitation index at inception versus advance coefficient. 
At the end of a day's testing, the water was drained from the 
test section. The propeller was operated in air, tare 
loadings, atmospheric pressure and no-flow manometer height 
were once again recorded. All raw data recorded is contained 


in Appendix B. 
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Vs DATA REDUCTION 


Reduction of the data for obtaining propeller parameters 
was accomplished using a program written for a TI-59 
programmable hand calculator. This program performed the 
following functions: 

(1) Determined the changes in tare loadings, 
atmospheric pressure and no-flow manometer height between 
the beginning and end of a series of test runs. A linear 
interpolation, based upon run number of a series, was then 
used to determine the value of these parameters for each run. 

(2) Air and water temperatures were used to determine 
the vapor pressure of water at the two temperatures (Di. and 
pe) The vapor pressure of water at the room air temperature 
was used to correct the reading of the ‘mercury column which 
was used for indication of tunnel static pressure, since this 
reading (p) was actually static pressure minus the vapor 
pressure of water at room temperature. Tunnel water temperature 
was also used to determine its density, op, and kinematic 
viscosity, v. 

(3) Tabulated values for the conversion from manometer 
height to velocity for the range of values encountered ina 
given test were entered and stored. They were then used in a 
linear interpolation to determine flow speed. 

(4) For each value of static pressure where data was 
recorded during the test, the following calculations were 


made: 
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(a) Manometer height was entered, corrected 
for the no-flow zero and converted to flow speed, 
U., in feet per second. 

(b) Static pressure reading, p, was entered. 
To this value the vapor pressure of water at room 
temperature, Pia? was added to give the true static 


pressure at the tunnel axis: 


Pstat = 2 Pyva 


(c) Propeller RPM was entered and converted to 
revolutions per second, n. 

(d) Propeller thrust reading was entered and 
corrected for the tare loading. A correction for 
the change in thrust from the pressure differential 
between tunnel static pressure and atmospheric 
pressure acting across the 1.317 in diameter 
propeller shaft was then applied to give the actual 
Encuste, ©. The Ehrust coefficient 

As 
on*p* 
where D is the propeller diameter, was then calculated. 

(e) This value of Ky was used to determine 
the advance coefficient, J, from the open water test 
results provided for this propeller, This value 
was entered and stored for later use. 

(£) For the first pressure increment for each 


test, the measured torgue was then entered, where 





a= 


it was corrected for the tare loading to give 


Measured torque, Q, and a torque coefficient, 


was calculated. The open water test results 

were entered with this value of Ky to verify that 
the J obtained in the previous step was reasonable. 
The value of J obtained from the thrust identity 
was used for the cavitation inception curve. In 


addition, at this step, the 0.7 radius Reynolds 


number, 


L. 
C x (V_°+(0.77nD) *) * 


6.7 
Vv 


: where We is the average inflow velocity seen at the 


propeller (calculated from J) andc was the 


0.7 
blade chord at the 0.7 radius, was calculated. 


(g) Finally, the cavitation number, 


7 Pstat ~ Pyw 
4US 


Once the data was reduced, and a value of J and o could 
be asSigned to each data run, it was necessary to determine 
which run represented the cavitation index at inception, Oss 
for each value of J. The criteria used to define inception 
are as follows: 

(1) Visual observation - Hub vortex cavitation had a 


different criterion from all other types. For it, the 
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appearance of a trail of bubbles from the vicinity of the hub 
was used as the criterion. For other types of cavitation, the 
criterion was to have a steady occurrence of that type of 
cavitation on more than one blade. Steady occurrence meant 
that the cavitation was present on each revolution of the 
propeller at one location, but not necessarily throughout 
the entire revolution. 

(2) One third octave band level - the arbitrary db 
level displayed on the measuring amplifier meter was plotted 
against decreasing cavitation index as shown in figure 8. 

The value of o which corresponded to the curve after the "knee" 
being 3 db above the extension of the curve before the knee 
was taken as O;. 

(3) Spectral analysis, complete spectrum - an increase 
of 3 db from the level at atmospheric pressure across the 
40-50 kHz portion of the spectrum was taken as the criteria. 
The three spectra in figures 9, 10 and 11 correspond to a 
fixed J, 0.51, and three different values of o: that for 
atmospheric pressure, o. for acoustically determined inception 
and Oo. for acoustically determined inception, respectively. 
The differences between these spectra are very slight and 
tend to make the determination of inception difficult and 
rather arbitrary. For this reason, this method was abandoned 
in favor of using the demodulated analysis. 

(4) Spectral analysis, demodulated signal - it was 
assumed that the presence of a sharp peak ("line") at shaft 


rate frequency indicated one blade cavitating, and that a 
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line at blade rate (number of blades times shaft rate) 
indicated all blades were cavitating (despite the possibility 
of the blade rate line merely being a harmonic of the shaft 
rate line). The inception criteria for this analysis was 
first taken to be the presence of lines at both shaft rate 
and blade rate which were at least 3 db above the general 
trend of the noise. A subsequent and less stringent criteria 
finally adopted was to require the presence of a line 3 db 
above the noise at blade rate frequency, with other lines 
present at shaft rate spacing to verify that cavitation was 
causing the line. Typical demodulated spectra are shown in 
figures 12 and 13. 

For each test, the value of Oo; obtained was plotted 
against J to produce the cavitation inception curve. The 
two different inception criteria for the demodulated signal 
analysis were plotted on the same graph, but with 


contrasting symbols and curves. 
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V. RESULTS AND DISCUSSION 


Curves of o, versus J for the visual and acoustic 
cavitation inception determinations are presented in 
figures 14 through 18. The comparison of acoustic and 
visual determinations for the one-third octave acoustic 
level measurements (figures 14 and 15) show a reasonably 
good agreement between the acoustically determined curve 
and the portion of the visually determined curve corresponding 
to the back bubble cavitation. However, tip vortex and 
leading edge face cavitation do not appear to be detected 
acoustically with the setup used. 

Comparing the acoustic and visual results for the 
demodulated analysis of the acoustic signal figures 16, 17 
and 18 shows a much better agreement for all types of 
Cavitation, except for hub vortex cavitation. The results 
Of the demodulated analysis of the 50-63 kHz band (figure 17) 
agrees almost exactly with the visually determined results 
for both of the inception criteria used with the acoustic 
analysis. The demodulated analysis of the acoustic signal 
above 20 kHz (figure 18) shows cavitation inception occurring 
at a higher value of o for tip vortex cavitation than for 
the visual results, and at about equivalent values of o for 
leading edge cavitation. The less stringent criteria for 
acoustically determined inception shows cavitation occurring 
at a higher value of o than the more stringent criteria. 

The two sets of results presented use different 


acoustic sensors. For the one-third octave level measurements, 
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Figure 15 


Acoustically Determined 
Cavitation Inception 
(Inception based upon one- 
third octave band level. 
Data point iS average over 
three bands - 30 to 60 KHz) 
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Figure 16 
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Figure 18 
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the accelerometer was used; but the demodulated analysis 
presented used the hydrophone, since it was felt that the 
higher usable frequency for the hydrophone was necessary. 
However, a test run at J=0.62 (run number 2 of 4 March) was 
performed to compare the acoustic information obtained from 
the two sensors. Examples of the demodulated spectra from 
the two sensors for a given o are shown in figures 19 
(hydrophone) and 20 (accelerometer). Except for the equipment 
gain adjustments needed to accomodate the different 
sensitivities of the sensors, the spectra are almost 
identical, indicating that either sensor was usable for an 
acoustic detection method. 

It was expected that the curves of oF versus J would 
show good agreement between the acoustic and visually ‘ 
determined cavitation inception, and these results confirm 
this. It was also expected that the acoustically determined 
inception would anticipate (occur at a higher value of o) 
the visually determined inception. This, in general, did not 
SCCur . 

The higher value of o at acoustically determined 
inception is based upon bubble size considerations. It was 
first assumed that the minimum bubble diameter which could 
be detected visually under the conditions of a propeller 
Cavitation test was 0.001 in. This size was then taken to be 
the maximum diameter (2R)) of a bubble in the calculation 


Shown by Strasberg (1977) for the total lifetime of the cavity, 


T= 2.7 Ry 1 OEE 
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Figure 19 
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where op and Po were taken to be representative of the 
conditions in a cavitation test, ‘p = 1.93 lb-sec’/ft’, 
P, = 400 mm Hg, or 1114 iy eto ender sehese conditions, 
T = 5.62 x 10 * sec, which corresponds to a frequency peak, 
4! in the cavitation noise spectrum of 17 kHz. With Po 
equal to atmospheric pressure, =. = 24.5 kHz. Since the 
frequency used for this analysis was much higher than these, 
it was felt that the detection of inception would occur at a 
higher value of o. 

Two possible explanations for the observed result not 
being in agreement with the prediction come to mind. The 
first is that the existence of scale effects (affecting the 
frequency scaling), due particularly to compressibility, 
surface tension and viscous effects,,were not ‘ttaken into 
account. If this caused the discrepancy, the use of high 
frequency acoustic information to anticipate visual inception 
determination would not be a workable scheme. However, if the 
expected acoustic signal was present, but was not detectable 
With the method or equipment used here, then anticipating the 
Visual inception determination is possible, so long as the 
appropriate changes are made. 

If the 0.001 in diameter bubble mentioned above is used 
With Strasberg*s non-dimensionalization for acoustic power, 
the ratio of the power output with the spectral peak at 56 kHz 


to that with the just barely visible bubble are given by, 


‘Oe 


acCOuStlIc a. V 17 


JO ag 


visual 7, a 
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or about 10 db. Thus, without considering the noise present 
or the increased absorption of the high frequency signal, 
10 db of gain are required to have an equal acoustic signal 
with the two conditions. When these other considerations 
are included, it becomes obvious that increasing the gain of 
the signal or decreasing the level of noise, or both, is 
needed. The dramatic change between the inception information 
obtainable with a one-third octave analysis and that 
obtainable with the demodulated analysis tends to verify this. 

By using an acoustic sensor which had some degree of 
directivity, either by using an array of hydrophones or by 
using some sort of reflector, an increase in the signal 
to noise ratio could be expected. Problems encountered ° 
with the instrumentation used could be corrected: 

(1) There was no account taken of the changes in 
absorption that occurs as air bubbles grow when the pressure 
is reduced in the test section. The use of a calibrated 
reference signal in the frequency range of interest would 
enable correction of the acoustic signal levels for absorption. 

(2) The combination of spectrum analyzer and X-Y 
plotter used required about three minutes to produce a paper 
copy of the demodulated spectrum, and the concurrent acoustic 
and visual cavitation inception determination extended this 
time span to the range of four to five minutes. Thus, for 
Gach value-of static pressure for a given J, about four 
Minutes were required, and the time required to produce each 


data point on the o, versus J curve was about thirty minutes. 
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Because of time constraints on the availability of the test 
facilities and equipment, and the time required for each data 
point, two factors added to the inaccuracies in the results. 
First, each data point represented only one test at that 
value of J. Second, the steps in tunnel static pressure used 
were on the order of 25 mm Hg. This represents an error of 
2.units.of -cavitation index at a tunnel flow speed of 6 feet 
per second or 0.3 units at 16 feet per second. 

A more rapid analysis of the demodulated spectrum would 
make the method less time-consuming. 

(3) During the time required for the spectrum analysis 
and averaging, tunnel flow speed and static pressure and 
propeller RPM would tend to drift on the order of one to five 
percent., The Benue cine effects of these changes would also 
affect the accuracy of the analysis by causing variations 
(although slight) in the frequency of interest and by 


affecting the values of o and J for the test run. 





ise 
VI. CONCLUSIONS 


Strasberg (1977) points out that "it is not possible 
to estimate the inception cavitation number of the prototype 
from model measurements without using empirically or 
theoretically determined scale factors." However, the results 
here show that it is possible to determine the cavitation 
inception performance of a model propeller by acoustic means 
at least as accurately as by visual means, as long as an 
adequate system for detecting the noise from all types of 
cavitation was available. And although the acoustically 
determined inception would require the same scale factors 
mentioned by Strasberg to predict full scale inception, the 
use of an acoustic inception determination technique for 
model tests does have advantages. 

First, where full scale inception measurements are 
made acoustically, an acoustic measurement technique for 
the model would eliminate any scale effect that would occur 
between model and full scale measurements caused by visual 
observation on the model and acoustic determination on the 
full scale propeller. Although the results indicate that for 
the propeller tested here this scale effect would be small, 
the test of a different propeller, with a different length 
scale might show a greater difference between acoustic and 
visual results. 

Second, displaying the spectrum of the demodulated 
Cavitation noise signal gives a more definitive criteria for 


inception than visual methods, as expected. 
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Appendix A. 
Details of the Wake Screen Design 


The method developed by McCarthy (1963) and adapted into 

a FORTRAN program by Rose (1969) attempts to determine the 
value of the non-dimensional resistance coefficient for a 
mrid, 

K = a i) 

oes 

where Ap is the local static pressure drop across the screen 
grid and Wo is the local velocity normal to the grid. 
McCarthy points ame Gaede an empirical estimate for K for a 


given screen laid over a support screen is given by 


S 
het 
Oe fs: Paar 2 (2) 
where: s, the solidity ratio for the screen = 


MD (2-MD) 

M is the number of wires per inch in the mesh 

D is the diameter of the wire in inches 

K ig the resistance coefficient for the support 

screen 

The program written by Rose requires that the test section 
area be subdivided into smaller areas, Ass with a flow 
Velocity, Vas associated with each area, which is the average 
Velocity for the subdivision. The overall velocity average 


is then calculated, 





S15 





V = role (3) 


Then, for each area, the resistance coefficient, Key 1s 
calculated as follows: 

(1) The integration constant, ie from McCarthy's 
solution is determined from the area having the maximum 
average velocity, V , which is assumed to have the resistance 


MAX 


coefficient, Kor for the support screen only: 





eS 2 \1/3 
y= = eee _Ymax_ -~1{+y_- a ie) 
O S 
N X tl avg 6 (K_+1) 
where y= Clk \% 
S S 
N = 1.02 : 


(2) For each area, A. with its associated Vir K. is 


determined by solving the following equation for Ky: 





Ua 1 Xs +1 1/3 
ee ee | N a. > (5) 
“ae 6 (K, +1) (2+K. +x, ) 
h = (1+K 2 
where es ; 
N =1.02 


(3) Once Kis determined for each area, the last step 
is to account for the deflections of the streamlines by the 
Changes in velocity caused by the wake screen. This is done 
by an iterative procedure which adjusts the areas, A, Ehgieat ak 
the volume flow rate at the screen is equal to the volume flow 


Yate at the propeller (this assumes that the screen is within 
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one tunnel diameter of the propeller, and it uses an empirical 
constant, a4). The result is a correction to the actual screen 
area for each subdivision. 

For this experiment, this program was adapted for use 
with a TI-59 programmable hand calculator. The test section 
was divided into four areas, each with its associated velocity, 
aS shown in Table A-l. These areas and velocities were based 
upon the desired circumferential mean wake information provided 
for the velocities over the propeller disk, and an estimate for 
the average velocity outside that area. For each area, Ks and 
the area correction, AA.» were determined from the program. 

At this point Rose's method and the one used here 
become different. When a large number of screens with different 
Meshes mad wire diameters are available, having a screen 
available with the proper resistance coefficient enables the 
Final wake screen to be assembled by piecing together the 
correct screens on the support screen. But here, the desired 
screens were not all available, so an alternative had to be 
developed. 

The alternative was based on an interpretation of 
information given in Pope and Harper (1966) on turburlence 
generation by screens in wind tunnels. This text indicates 
that the cumulative effect of several layers of screens was 
additive. This seemed to be supported by the empirical 
formula for K in equation (2), where the effect of the epost 


screen and the wake producing screen are added. It was 
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Table A-l 


Wake Characteristics 





Area Range of A Vy A r/R 
- . K 
Designation r/R i ot i CO corr req 
1 eo =- 1,0 308.39 .97 317.86 o =—-,946 0.4612 
2 1.0-0.7 46.72 .817 44.51 .946-.641 0.9462 
3 ree ora 299 V656 1827S .641=—.454 1.756 
4 -5-.23 omc ahora. Ole. 454—.23 . 2.531 
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felt that, even though Pope and Harper note that the effect of 
the screens was additive only if they did not touch, assuming 
that the effect was additive would be satisfactory for a first 
approximation. 

Thus, the desired initial K values were obtained by 
using several layers of two different screens. The 
characteristics of these screens are listed in Table A-2, and 
the actual Ks for the screens used is listed in Table A-l. 

The was screen was then assembled, with the pieces of screen 
wired onto the support screen with pieces of 0.020 in stainless 
steel wire. A wake survey conducted with the Laser Doppler 
Anemometer was performed along the diagonal line shown in 
figure 3 previously. 

The determination of velocities in the wake was done 
with the propeller removed from the shaft. The longitudinal 
Flow speeds were determined in the plane 2.5 in downstream of 
the propeller blade leading edges. The results of this and all 
Other wake surveys performed are contained in Appendix B.1. 
The results for this initial screen are plotted in figure A-l. 

From the plot of non-dimensional velocity versus non- 
dimensional radius, r/R, it was possible to determine the 
average velocities actually obtained for each of the 
Subdivided areas used, since the wake was to be axisymmetric. 
From the actual values of Vir equations (4) and (5) were then 
used to determine the actual K obtained from the screen used 


for each area, shown in Table A-3. At this point it was 
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Table A-2 


Screen Characteristics 


Solvdity 

Wire Diameter, D ratio, s 
0.020 in 2944 
O02 an 28853 


Table A-3 


Kan 
= 


-4612 


a 956 


Comparison of Installed and Measured K Values 


Area Kreg 
1 0.4612 
2 0.9462 
3 1.756 
4 yee 3 1 


ee 


AK ast Kneas AKneas 
0.4612 - 0.4612 - 
0.9224 0.4612 pa Ce es. O27 20) 
1.7179 0.7954 2.4197 1.2384 
2.5134 0.7954 4.106 1.6863 


AK ratio 


JE Baie) 0, 
JER SS 7 
Za 0 
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assumed that the calculated drag associated with the support 
screen was correct, and the remaining screen layers could 
be adjusted to develop the desired wake. 

Still assuming that the effect of multiple layers of 
screen material was additive, the measured increase in K EO 
each fe rement of screen material was determined (AK, in Table 
A-3). In each case, the ratio of measured AK, to installed 
AK was determined. The results of these calculations were 
interpreted as follows: 

(1) In area 2, the same material as that found in 1 was 

added, and the increase in K was 1.56 times greater than 
that macted from the simple addition of resistance 
coefficients. 

(2) Comparing measured results for areas 3 and 4, the 

increase from adding another layer of the same material 


1.6863 _ 1.36) times greater than the expected 


was 1.36 ( T.2354 = 


result. 

(3) The apparent effect of adding the screen of area 3 

onto that of area 2 was 1.56 greater than expected. 

(4) Thus, an average increase in the resistance 

coefficient from adding layers of screen material was 

Felt to be on the order of 1.5 times the expected result. 
In terms of the screen material available, this meant that the 
8x8x.020" screen material would be used for all screen layers, 
with each added layer having 1.5 times the resistance 


coefficient of the single layer. Table A-4 shows the 
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O- measured 





Wake Survey - Initial Wake Screen 


Figure A-l 





Wake Survey - Final Wake Screen 


Figure A-2 





-65- 


calculated eeiaigied 


compared to the required K values, where 


values for the four wake screen areas 


This screen was assembled as noted previously, and another 
wake survey was Made at the same points as before. 

The results of this survey are shown graphically in 
figure A-2. The low velocity seen at the outer radii 
correspond to the use of a screen giving higher than the 
required value of K in area 2. Lacking a screen with a 
small enough K to correct this discrepancy, this wake was 


considered to be acceptable. 


Table A-4 


Final Screen Resistance Coefficients 


Area Kreg K added KR wale 
1 0.4612 O24 6.22 
2 0.9462 04 Geb2 148 
3 7 6 0.4612 1.836 
4 2a 0.4612 Zee Ss 
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APPENDIX B 
Raw Data 


The following pages of this section contain the propeller 
Operating conditions data sheets, the graphs of sound level 
versus cavitation index for one-third octave level measurements, 
and the X-Y plotter outputs for the demodulated spectra and 
for the analysis of the complete spectrum. Since, in some 
cases, the recorded data may not be completely clear as to 
what is being presented, some clarifying explanations are 
presented here. 

Wake survey data - For tables B-l and B-2,. the velocity 
used to obtain (l-w) is the Veavg calculated from 
the 6,7,and 8 inch radial positions. For table 
B-3, the velocity used is the velocity from the 
Manometer height, 12.0 fps. 

One-third octave band level - The conversion factors for 
T to thrust and Q to torque for this series of tests 
and for the data of sections B.5, B.6 and B.7 are on 
the data sheet for run number 1. In these and all 
subsequent data sheets the first number in the 
"GAIN CHANGE" column refers to the code number for 
the amplifier shown in the upper section, the second 
number gives the setting to which the selector was 
changed. In the "REMARKS" column for this section 


only are three numbers which refer to the displayed 
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level on the measuring amplifier meter for the 
31.5-40 kHz band, the 40-590 kHz band, and the 50-63 
kHz band,respectively. For all sections of the 
appendix this column also contains the visual 
inception determinations as follows: any visual 
inception determination will have a * , together 
With a LV £Or tip vortex cavitation, HV for hub 
vortex cavitation, LEPF for leading edge pressure 
face cavitation, LESF for leading edge suction face 
cavitation, BB or BACK for back bubble or back 
cavitation. The points plotted on the plots of db 
versus oO are the arbitrary level from the measuring 
amp meter. The top plot is the aS =40 kHz level, 
the middle is the 40-50 kHz level, and the bottom 
is the 50-63 kHz level. 

Demodulated Analysis - The conversion factors for thrust 
and torgue for these two sections are on the first 
data sheet for the 20 kHz high pass signal. On 
these and all plots from the X-Y plotter, 10 db is 
represented by 20 of the smallest divisions (2 cm 
total). For the plots in these two sections, the 
first plot made was the one at the bottom of the 
page. For each subsequent plot, the zero setting 
for the X-Y plotter was placed 2 cm higher (10 db). 
Unless noted on the associated data sheet, no gain 


adjustments were made. On these and all subsequent 
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plots the number typed in beside a particular curve 
is the value of the cavitation index associated with 


that plot. 
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Table B-lL 


Wake Survey Data - 


Initial Wake Screen 


Scale 

Hor. Vert. 
11.74 91.60 jreislly 
F300. 39.390 Oo 7, 
14.43 88.01 Bec Pe ler 
ES. 24.66.93 Vie 200) 5 
Poo oo. 96 Sos? , 
16.69 84.99 yO22; 
17.42 84.02 .546, 
Poets oo Oo -443, 
18.87 82.08 Scio ar 
AMSHA SIS) roel ee E 200; 
20.33 80.14 5 FOE FE 
20.84 79.46 28+ 
Zoe s Fo.00 SIM er 
24.69 74.32 50; 
25242 73.235 73055 
26505 J2<358 -408, 
26.87 71.41 ~443, 
27.60 70.44 -540, 
28.33 69.47 2O24, 
ASS) 610 ey (ESS SiO Gas, 
29.78 67.53 »/159, 
30.59 66.45 .844, 
31.94 64.66 Sissel - 
De Te 5 2 
= 20.805 ft/sec 

836 V = 17.4 £t/sec 
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tT © 
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(l-w) 
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Sua 
Aiea: 
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- 436 
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Be) 
~354 
- 394 
-408 
~ 433 
-487 
5S) 3)10, 
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~/44 
ei) 
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Table B-l (Continued) 


Comparison With Prescribed Wake 


Abe aby Wake Screen 


r/R (Lew) avg WET oye Baeror 
iiss 5 903 -.03 
O59 5 ES) .846 -~.10 
0.8 ~ 745 - 784 =~ 05 
0.7 .649 .710 -.09 
Giro 529 7036 ~.17 
0.5 -461 273 -.19 
0.4 ~385 SICA: ~.27 
02.3 305 488 ~.25 


0.23 -374 -431 Se dns 
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Table B-2 


Wake Survey Data - 


Final Wake Screen 


Scale 
Hor. Vert. 
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Table B-2 (Continued) 


Comparison With Prescribed Wake 
Final Wake Screen 


(Unmodified upper support) 


r/R lew) (ED ys Error 
IEA 6 ~901 ~903 =-.003 
0.9 ~ 746 .846 -.118 
0.8 e732 -784 -.066 
O57 -659 = 7 Auo O72 
0.6 -593 O36 -.068 
05 .950 257-5 -.040 
0.4 ~491 ~528 -.070 
0.3 -475 -488 a 024! 
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Table B-3 


Wake Survey Data - 


Final Wake Screen 
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Table B-3 (Continued) 
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Scale 


Hor ® 


24.58 
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Table B-3 (Continued) 
Volts 
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Table B-3 (Continued) 


y Zz Scale Volts Speed (l-w) 
Hor. Vert. 

-4 4 £Se06 68. 22 ~629 13.09 i. 09 

3 70.76 ~495 1023.0 . 86 

2 737. S10 ~459 9.55 730 
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-1 80.92 soo. 13.54 Les 

-6.5 -6.5 14.10 94.89 tO Ou 2.54 1s 

-6 -6 Poo 95.02 woos 159 pots 

=5.5 =-5.5 16.01 92.35 .666 13.86 5 

—5 -5 16.96 91.08 FOO 7 VSe67 ee 
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Coordinates are in inches. (l-w) is based on 12 fps 
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